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Abstract 

We present the results of a complete leading-twist next-to-leading order 
(NLO) QCD analysis of the spacelike pion electromagnetic form factor at 
large momentum transfer Q. We have studied their dependence on the 
form of the pion distribution amplitude. For a given distribution ampli- 
tude, we have examined the sensitivity of the predictions to the choice 
of the renormalization and factorization scales. Theoretical uncertainty 
of the LO results related to the renormalization scale ambiguity has been 
significantly reduced by including the NLO corrections. Adopting the cri- 
teria according to which a NLO prediction is considered reliable if, both, 
the ratio of the NLO to LO contributions and the strong coupling con- 
stant are reasonably small, we find that reliable perturbative predictions 
for the pion electromagnetic form factor with all distribution amplitudes 
considered can already be made at a momentum transfer Q < 10 GeV, 
with corrections to the LO results being typically of the order of ~ 20%. 
To check our predictions and to discriminate between the distribution am- 
plitudes, it is necessary to obtain experimental data extending to higher 
values of Q. 
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I. INTRODUCTION 



Exclusive processes involving large momentum transfer are among the most interest- 
ing and challenging tests of quantum chromodynamics (QCD). 

The framework for analyzing such processes within the context of perturbative QCD 
(PQCD) has been developed by Brodsky and Lepage |J, Efremov and Radyushkin 0, 
and Duncan and Mueller || (see Ref. [|J for reviews). They have demonstrated, to 
all orders in perturbation theory, that exclusive amplitudes involving large momentum 
transfer factorize into a convolution of a process-independent and perturbatively incal- 
culable distribution amplitude, one for each hadron involved in the amplitude, with a 
process-dependent and perturbatively calculable hard-scattering amplitude. 

Within the framework developed in Refs. [0-0], leading-order (LO) predictions have 
been obtained for many exclusive processes. It is well known, however, that, unlike in 
QED, the LO predictions in PQCD do not have much predictive power, and that higher- 
order corrections are essential for many reasons. In general, they have a stabilizing effect 
reducing the dependence of the predictions on the schemes and scales. Therefore, to 
achieve a complete confrontation between theoretical predictions and experimental data, 
it is very important to know the size of radiative corrections to the LO predictions. The 
list of exclusive processes at large momentum transfer analyzed at next-to-leading order 
(NLO) is very short and includes only three processes: the pion electromagnetic form 
factor [p[-{T0f1, the pion transition form factor flO-Ii^], and photon-photon annihilation 



into two flavor- nonsinglet helicity-zero mesons, 77 — > MM (M = tt, K) (13 

In leading twist, the pion electromagnetic form factor (the simplest exclusive quan- 
tity) can be written as 

FAQ 2 )= f dx f 1 dy &{y^ 2 F )T H {x,y,Q 2 ^ 2 R ^ 2 F ) <2>(x,/4). (1.1) 
Jo Jo 

Here <P(x,fj, 2 F ) is the pion distribution amplitude, i.e., the probability amplitude for 
finding the valence qiq 2 Fock state in the initial pion with the constituents carrying the 
longitudinal momentum xP and (1 — x)P; Th (x , y , Q 2 , fi R , fi 2 F ) is the hard-scattering 
amplitude, i.e., the amplitude for a parallel q{q 2 pair of the total momentum P hit by 
a virtual photon 7* of momentum q to end up as a parallel q{q 2 pair of momentum 
P' = P + q; $*(|/,/^) is the amplitude for the final state q{q 2 to fuse back into a pion; 
Q 2 = —q 2 is the momentum transfer in the process and is supposed to be large; [i R is the 
renormalization (or coupling constant) scale and fxp is the factorization (or separation) 
scale at which soft and hard physics factorize. 

The hard-scattering amplitude T H can be calculated in perturbation theory and rep- 
resented as a series in the QCD running coupling constant as{fi R )- The function $ is 
intrinsically nonperturbative, but its evolution can be calculated perturbatively. 

Although the PQCD approach of Refs. [|l]-|3j undoubtedly represents an adequate 
and efficient tool for analyzing exclusive processes at very large momentum transfer, its 
applicability to these processes at experimentally accessible momentum transfer has long 



been debated and attracted much attention. The concern has been raised |fl^fl5[1 that, 
even at very large momentum transfer, important contributions to these processes can 
arise from nonfactorizing end-point contributions of the distribution amplitudes with 
x ~ 1. It has been shown, however, that the incorporation of the Sudakov suppression 
effectively eliminates these soft contributions and that the PQCD approach to the pion 
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form factor begins to be self-consistent for a momentum transfer of about Q 2 > 4 GeV 2 
|T§ (see also Ref. |7|). 

To obtain the complete NLO prediction for the pion form factor requires calculating 
NLO corrections to both the hard-scattering amplitude and the evolution kernel for the 
pion distribution amplitude. 

The NLO predictions for the pion form factor obtained in Refs. [0-0] are incomplete in 
so far as only the NLO correction to the hard-scattering amplitude has been considered, 
whereas the corresponding NLO corrections to the evolution of the pion distribution am- 
plitude have been ignored. Apart from not being complete, the results of the calculations 
presented in Refs. [[|-0] do not agree with one another. 

Evolution of the distribution amplitude can be obtained by solving the differential- 
integral evolution equation using the moment method. In order to determine the NLO 
corrections to the evolution of the distribution amplitude, it is necessary to calculate two- 
loop corrections to the evolution kernel. These have been computed by different authors 
and the obtained results are in agreement fl8|] . Because of the complicated structure of 
these corrections, it is possible to obtain numerically only the first few moments of the 



evolution kernel fl9|| . Using these incomplete results, the first attempt to include the 
NLO corrections to the evolution of the distribution amplitude in the NLO analysis of the 
pion form factor was obtained in Ref. ||10|| . It has been found that the NLO corrections 
to the evolution of the pion distribution amplitude as well as to the pion form factor are 
tiny. 

Considerable progress has recently been made in understanding the NLO evolution of 
the pion distribution amplitude pD |. Using conformal constraints, the complete formal 



solution of the NLO evolution equation has been obtained. Based on this result, it has 
been found that, contrary to the estimates given in Ref. ||10|| , the NLO corrections to 
the evolution of the distribution amplitude are rather large. It has been concluded that 
because of the size of the discovered corrections, and their dependence upon the input 
distribution amplitude, the evolution of the distribution amplitude has to be included in 
the NLO analysis of exclusive processes at large momentum transfer. 

The purpose of this paper is to present a complete leading-twist NLO QCD analysis 
of the spacelike pion electromagnetic form factor at large momentum transfer. 

The plan of the paper is as follows. To check and verify the results obtained in Ref. 
in Sec. II we carefully calculate all one-loop diagrams contributing to the NLO 
hard-scattering amplitude for the pion form factor. We use the Feynman gauge, the di- 
mensional regularization method, and the modified minimal-subtraction (MS) scheme. 
Our results are in agreement with those obtained in Refs. [0] (modulo typographical 
errors listed in ||) and ]7). Making use of the method introduced in Ref. p0| , in Sec. Ill 
we determine NLO evolutional corrections to four available candidate pion distribution 
amplitudes. In Sec. IV we discuss several possible choices of the renormalization scale 
fin and the factorization scale /ip. In Sec. V we obtain complete NLO numerical pre- 
dictions for the pion form factor using the four candidate pion distribution amplitudes. 
For a given distribution amplitude we examine the sensitivity of the predictions on the 
renormalization and factorization scales, /ir and /if, respectively. We take Aj^g = 0.2 
GeV for the calculation presented here. Section VI is devoted to discussions and some 
concluding remarks. 
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FIG. 1. Feynman diagrams describing the (qiq 2 ) + 7* ~~ * (Q1Q2) transition amplitude in 
terms of which the hard-scattering amplitude for the pion form factor is obtained. 

II. NLO CORRECTION TO THE HARD-SCATTERING AMPLITUDE 

In this section we recalculate the NLO correction to the hard-scattering amplitude 
for the pion form factor. 

In leading twist, the hard-scattering amplitude is obtained by evaluating the (qiq 2 ) + 
7* — * (Q1Q2) amplitude, which is described by the Feynman diagrams in Fig. [l|, with 
massless valence quarks collinear with parent mesons. In this figure, uP and (1 — u)P 
{yP' and (1 — v)P') denote the longitudinal momenta of the pion constituents before 
the subtraction of collinear singularities. In this evaluation, terms of order m 2 /Q 2 are 
not included and, since the constituents are constrained to be collinear, terms of order 
k 2 ± /Q 2 (k± is the average transverse momentum in the meson) are not taken into account 
either. By projecting the q{q 2 pair into a color-singlet pseudoscalar state the amplitude 
corresponding to any of the diagrams in Fig. [I] can be written in terms of a trace of a 
fermion loop. 

By definition, the hard-scattering amplitude T# is free of collinear singularities and 
has a well-defined expansion in as(/i|.) of the form 



T H (x,y,Q 2 ,fi 2 R ,fi%) = a s (^ 2 R )T^\x,y,Q 2 ) 

' , (2.1) 



x 



i + *lr«(x,^ 2 fl /QVW) + 



7T 



where 



and 



4tt 

asi ^ = MAiKl CD y (2 - 2) 



A = U-|n/, (2.3) 

with rif being the effective number of quark flavors and Aqcd is the fundamental QCD 
parameter. 

In the LO approximation (Born approximation) there are only four Feynman dia- 
grams contributing to the (qiq 2 ) +7* - ► {Q1Q2) transition amplitude. They are shown in 
Fig. 0. Evaluating these diagrams, one finds that the LO hard-scattering amplitude is 
given by 

Tg\*,y,Q 2 ) = A , n2h 16 : n y (2-4) 

3Q 2 (1 -x)(l -y) 
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FIG. 2. Lowest-order Feynman diagrams contributing to the (q^) + 7* (Q1Q2) amplitude. 

To obtain this result, it is necessary to evaluate only one of the four diagrams. Namely, 
knowing the contribution of diagram A, the contribution of diagrams B, C, and D can 
be obtained by making use of the isospin symmetry and time-reversal symmetry. In the 
leading- twist approximation, the isospin symmetry is exact. A consequence of this is 
that the contributions of diagrams A and B are related to the contributions of C and 
D, respectively. On the other hand, diagrams A and C are by time-reversal symmetry 
related to diagrams B and D, respectively. 

At NLO there are altogether 62 one-loop Feynman diagrams contributing to the 
{Q1Q2) + 7* ~^ {Q1Q2) amplitude. They can be generated by inserting an internal gluon 
line into the leading-order diagrams of Fig. 2. Use of the above mentioned symmetries 
(isospin and time-reversal) cuts the number of independent one-loop diagrams to be 
evaluated from 62 to 17. They are all generated from the LO diagram A. We use the 
notation where Aij is the diagram obtained from diagram A by inserting the gluon line 
connecting the lines % and j, where i,j = 1, 2, • • • , 6. They are shown in Fig. 3 with the 
exception of A33, A55, and A66, which give obviously the same contribution as All. 
These diagrams contain ultraviolet (UV) singularities and, owing to the fact that initial- 
and final-state quarks are massless and onshell, they also contain both infrared (IR) and 
collinear singularities. We use dimensional regularization in D = 4 — 2e dimensions to 
regularize all three types of singularities, distinguishing the poles 1/e by the subscripts 
UV and IR (D = 4 — 2euv = 4 + 2em) . Soft singularity is always accompanied by two 
collinear singularities and, consequently, when dimensionally regularized, leads to the 
double pole l/e 2 R . 

Before the subtraction of divergences has been performed, the pion form factor con- 
volution formula reads 

FAQ 2 ) = C du C dv &(v) A(u, v, il 2 vv IQ\ tf R /Q 2 ) *(«) , (2.5) 
JO Jo 

where A(u, v, [x^jyjQ 2 ■, iAr/Q 2 ) denotes the amplitude for the (qiq 2 ) +7* — ► (Q1Q2) quark 
subprocess. This amplitude has the expansion of the form 

A(u, v, ti v /Q\ tf R /Q 2 ) = a s A<°> («, v, Q 2 ) + — A (1) (w, „, ^ UV /Q 2 , tf R /Q 2 ) + ... , 

71 

(2.6) 

where 
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FIG. 3. Distinct one-loop Feynman diagrams contributing to the {qiq?) + 7* 
amplitude. The total number of diagrams is 62. 



(9192) 
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AW(u,v J Q 2 ) = (l-e)TH\u J v,Q 2 ). 



(2.7) 



The contributions of the individual diagrams to A^\u, v, jJqj V /Q 2 i V 2 r/Q 2 )i w hh the 
overall factor T^\u,v,Q 2 ) of Eq. ( P^ ) extracted, are listed in Table |. We have used 
the following notation: 



u 

Vuv 

VlR 
VlR 



1 — u, v = 1 — V , 

i , Q 2 

— 7 — in 



1 



4:7Tfx 2 uv e uv 
Q 2 l 



+ In 



Vuv 



+ 7 + In =- 

f//,< ^VlR e IR 

Q 2 



-In 



Q 2 

iAr 
Q 2 ' 



e //J e IR 
+7 In 
1 



i2 



7T 

IF 



+ -ln 2 







- 7 + ln(47r) , 
— — = — — "I - 7 — ln(47r) . 

ClR ClR 

The function H(w, t>) appearing in Table |I| is given by the expression 

1 



(2.1 



H(u,v) 



u — V 



Li 2 ( - ) + Li 2 ( - 

u J \v 



. . I uv 
+ Li 2 — 



-Li 2 ( — J - Li 2 ( — J - Li 2 ( — 



(2.9) 



where Li 2 (w) is the Spence function defined as 

Ll 2 M = - / ; 



(2.10) 



Adding up the contributions of all diagrams, we find that the NLO contribution to 
the (qiq?) + 7* — * (<?i<? 2 ) amplitude, given in ( |2.6|) , can be written in the form 



^ l \u,v^ 2 uv lQ\tiRlQ 7 
= T§\u,v,Q 2 ) ' 

+fuv(u,v,fi 2 uv /Q 2 ) + f IR (u,v,n 2 IR /Q 2 ) + fc(u,v) 



CuvJ— + C IR (u,v) — 
£uv £ir 



where 



Cuv 
C m (u,v) 



[3 + ln(wU)] , 



(2.11) 

(2.12a) 
(2.12b) 



and 
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TABLE I. Contributions to A^fu, v, H 2 R /Q 2 , ^ 2 F /Q 2 ) (denned by (jH)) of Feynman dia- 
grams shown in Fig. || . 



All 

A22 
A44 
A23 
A12 
A56 

A14 
A45 
A15 

A16 

A34 
A13 



A36 
A35 



1 

12 
1 
6 



(VUV + VI R) 
(VUV -Inu) 



> \ 16 4 

;n f j (rjuv ~ ln(HU)) + — - -n f 



If 1+v . _ 

■x Vuv H In v 

o V v 



u 



u \ u 

r/uv + %VlR ( 1 "I l 11 ^ ) — ^ Inu + Inv H In' 4 u + 2 — In u \nv 

u J u u u 



1 

48 

■— [r/[/v + 2r//^(l — ln(uu)) — 2?7/r — 5 + ln(uu ) — In 2 u — ln 2 U — 2 Inu lnv 
48 



3 

16 



3r] UV + 2r] IR 1 + 



Inu 



, u In 2 u 
+ 1 + In - + 

v u 



+ 2- 



lnu Inu 



Y£ [3r/[/V + 4f ?/K - 1 + m ( u « )] 



15 

3 u 

1 

24 
1 



In u + In u + - In 2 u + In u In u 

/ Inu _ , _\ lnu _ lln 2 u 2 

??J-R + Vi R 1 "I r 2 m u + in v H h 2 in u + m u H h in u 

\ u I u 2 u 



+ -ln 2 u + 



Inn Inu 



+ 2 in u In u 



1 



1 



/ In u In u \ _ u + v _ _ ■'-,■■»_ 1 . ■■> _ 
? ?7"H 2 H 1 — 2 + m u + m u H in u In u H In u H In u 

V u v J uv 2u 2v 

H in u H — in ?tH — in v + in u in u H in u in u 



2u 



2(u-u) 2 



2 2 u 

4u 2 uH(u, u ) + (4u - 5u 2 + u 2 )(lnu + lnu) 



2 , JJv^U , ,o„. ,..2 , ..... , ...2,. „.3>^ 



+(6u - 5u - 2u + u )u h (2u - 4u +2v - 4uu + 5u v - v 

u ' v 



1 
3 
1 

24 



u 



U + V 



r/j R — Inu -\ In u I H Inu -\ In u H In u -\ In v -\ '■ — In u Inv 

\u v J u v 2u 2v uv 

( 1+u _ 1 + u _\ 1 + u _ 

^ViR + ViR 2 H In u H Inv — 2(1 — u — v)tl{u, v) H Inu 

\ u v J u 



+ 



1 + _ . 1+U 2- , Hi!, 2 - , U+V _ _ 

1 1 in n in u 



Inu + 



2?7 



ln z u + 



2u 



ln^U + 



8 



a 



uv 



fuv(u,v,fi 2 uv /Q 2 ) 

f IR (u,V,tf R /Q 2 ) = C IR {H. r 

fc(u,v) 



ln(-ut> ) + In 



Q 2 



1 



-ln(uv) 



In 



Q 2 



(2.13a) 



(2.13b) 



— -10 

12 1 



20 ln(ii v ) + In m In v 

+ In U In t> — In u In v — In U In v 
+(l-u- v)H(u, v) + R(u, v)\ . 

The function R(w, i>) is denned as 

1 



(2.13c) 



R(w, u) 



(m — v) 2 

+ 



[(2w 



v)(lnu + In v) 



2m; 2 - 2v 2 + Wuv -2v- Au 2 



+ (-2vu 2 - 2u 2 + 10uu 
— (w 2 + mm 2 )H(m, 17 ) . 



2m - 4v 2 



In v 
i 

f 

lnU 
) 

u 



(2.14) 



The contributions of the individual diagrams listed in Table | are in agreement with 
those obtained in Ref. 0] (up to some typographical errors listed in 0). Correspondence 
between our results and those of Ref. is established by multiplying the latter by 
{P + P'Y/Q 2 . 

It is easily seen from Eqs. Q2.11 - 2T4]) that in summing up the contributions of 
all diagrams to the exclusive (gig 2 ) + 7* — > (Q1Q2) amplitude, the originally present 
soft singularities (double 1/e 2 ^ poles) cancel out, as required. The final result for the 
amplitude contains 1/euv an d only simple l/e IR poles. 

We now proceed to treat these poles. When doing this, one has to exercise some 
care, since the subtractions have to be performed in such a way that the universality 
of the distribution amplitude as well as the universality of the coupling constant are 
simultaneously preserved. 

The IR poles in ( |2.11|) , related to collinear singularities, are such that they can be 
absorbed into the pion distribution amplitude. The universality of the distribution ampli- 
tude requires that the poles 1/em should be absorbed by some universal renormalization 
factor H8|,|2ll1 . The regularization-dependent terms represent "soft" effects and therefore 
we absorb them into the distribution amplitude along with the singularities. A crucial 
observation is that the structure of the collinearly divergent terms in Q2.ll ) is such that 
one can write 



Vt(v, y) <g> T$\x, y, Q 2 ) ® S(x - u) + 6{v - y) ® Tg>(x, y, Q 2 ) ® Vx{x, it) 
= C IR (u,v)T { x\u,v,Q 2 ) , 



(0), 



(2.15) 



where the function V\(x,u) is the one-loop evolution kernel for the pion distribution 
amplitude (see the corresponding expression in, for example, or f22||), while <S> denotes 
the convolution symbol defined as 



A{z)®B{z) 



dzA{z)B(z) 



(2.16) 
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F 



(l-u)P (l-x)P (l-y)F (l-v)P' 

FIG. 4. Pictorial representation of the procedure for absorbing the collinear divergences 
into the pion distribution amplitude fj.'p). 



Based on Eq. (|2.15|) , the NLO expression for the (qiq 2 ) + 7* — > (Q1Q2) amplitude, given 
by Q2.6| ), (|2.11[ - pm ), can be written in a factorized form 



a s A<°> (u, v, Q 2 - e IR ) + 2* A«(t*, v, tf R /Q 2 ; l/e IR ) 

7T 



6(v-y) + ^V 1 (v,y) (± + \n^f- 

* \ e IR VlR, 



_ as, 
asT^\x,y,Q 2 ) 

L 71 

S( x - u ) + —V x {x, u) { J- + In ^f- 
n Vm Kir, 



[l-e) + ^f£\x,y,n 2 F /Q 2 ;e IR ) 

7T 



(2.17) 



Formula ( p. 17] ) shows the essence of the procedure for the separation of collinear diver- 
gences, which is schematically represented in Fig. f|. The procedure consists in summing 
the effects of collinear divergences contained in the shaded rectangular region in Fig. |], 
as a result of which the outside quark (antiquark) with momentum uP {uP) ends up on 
the inside quark (antiquark) of momentum xP (x P). 

As for the UV poles, we renormalize them using the modified minimal-subtraction 
(MS) scheme. This is carried out by stating that as appearing in (|2.6|) is the bare 
unrenormalized coupling related to the renormalized physical coupling o.-jjg{jji 2 R ) by 



"ms(/4) Po 



7T 



4 \iuv 



In 



(2.18) 



It is important to note that the e — > limit can be taken only after the separation 
of collinear divergences and UV renormalization have been performed (taking the e — > 
limit in A^ ) before the subtraction of IR and UV singularities is equivalent to choosing a 
factorization scheme that does not respect the universality of the distribution amplitude 
(see ||) and the renormalization scheme in which the running coupling constant is not 
universal). As a result, we obtain the following expression for the NLO hard-scattering 
amplitude for the pion form factor: 



T H (x, y, Q 2 , /4, H 2 F ) = «ms(/4) t h\ x , y, Q 



1 + 



«Mg(/4) 
7T 



T£\x,y,ti R /Q 2 ,» 2 F /Q 2 



(2.19) 



where 
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and 



T£\x,y,n 2 R /Q 2 ,y? F /Q 2 ) 

= fuv(x } y } ^ 2 R /Q 2 ) + f IR (x,y } fi 2 F /Q 2 ) + f c (x,y) 



(2.20) 



fuv(x,y,fx 2 R /Q 2 ) = f uv (x,y,fj, 2 R /Q 2 ) + C UV 



- ln(x y ) + In 



fm(x,y,n F /Q 2 ) = fm(x,y,fi 2 F /Q 2 
= ~[3 + \n(xy)} 



^ln(xy) -ln^ 



fc(x, y) = fc{x, y) - C IR {x, y) 

= — [—34 + 12 ln(x y ) + In x In y 

+ In x In y — In a; In y — In x In y 

+(1 - x - y)R(x, y) + R(x, y)] . 



(2.21a) 



(2.21b) 



(2.21c) 



A few comments on the previously performed calculations are in order. Our final 
expression for the hard-scattering amplitude Th is in agreement with Ref. 0. The 
calculation of || contains a few typographical errors in the diagram-by- diagram listing 
as well as in the final expression (which are all listed in 0). Apart from this, the 
results from Ref. U agree with ours. In Refs. |],[| the contribution of the diagrams 
with propagator corrections on external quark lines were not taken into account. Also, 
in Ref. the subtraction of collinear singularities was performed in a way that is not 
consistent with the universality of the distribution amplitude (see the discussion in Ref. 
H). Finally, a thorough analysis of the results of Refs. |5|-[7j was performed in Ref. 0, 
but in obtaining the final expression for the hard-scattering amplitude Th, collinear and 
UV divergences were subtracted in a way that is not consistent with the universality of 
both the distribution amplitude and the running coupling constant. 



III. EVOLUTIONAL CORRECTIONS TO THE PION DISTRIBUTION 

AMPLITUDE 

The pion distribution amplitude ^>(x,fi F ), controlling exclusive pion processes at 
large momentum transfer, is the basic valence wave function of the pion. Its form is not 
yet accurately known. It has been shown, however, that the leptonic decay tt + — > /U + z/^ 
imposes on $(x, fi F ) a constraint of the form 

Given the form of 3>(x, fi F ), this relation normalizes it for any [a f . In ( |3.1|) , / w = 0.131 
GeV is the pion decay constant and nc{= 3) is the number of QCD colors. 

Instead of using <P(x,fi 2 F ) satisfying (|3TT|), one usually introduces the distribution 
amplitude (f>(x, ji 2 F ) normalized to unity, 
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and related to <&(x, fi F ) by 



1 

dx <p(x, n F ) = 1 , (3.2) 

o 



$(x, up) = 2 Jl— (j)(x, iip) . (3.3) 



Although intrinsically nonperturbative, the pion distribution amplitude (f)(x, up) sat- 
isfies an evolution equation of the form 



fM F -^-(j)(x, fXp) = duV(x,u,a s (iip)) <t>{u,n F ) , (3.4) 

in which the evolution kernel is calculable in perturbation theory: 
V(x, u, asifip)) 

a s {fJ, F ) t r / s . ( a s (fi] 



2 \ 1 (x,u)+(^^-) V 2 (x,u) + 



IX y 7T 

(3.5) 

and has been computed in the one- and two-loop approximations using dimensional 
regularization and the MS scheme. 

If the distribution amplitude <f>(x, //q) can be calculated at an initial momentum scale 
/Xq using QCD sum rules |23[] or lattice gauge theory |24|, then the differential-integral 
evolution equation (|3.4j ) can be integrated using the moment method to give (p(x,fip) at 
any momentum scale fip > /ig- 

Because of the complicated structure of the two- loop contribution to the evolutional 
kernel V-^x, y), only the first few moments of the evolutional kernel have been computed 
numerically. 

Recently, based on the conformal spin expansion and the conformal consistency re- 
lation, the analytical result for the evolution of the flavor-nonsinglet meson distribution 
amplitude has been determined f20| . 

To the NLO approximation, this, rather complicated solution has the form 

<t>(x, /4) = <j> LO (x, + <f> NLO (x, lip) , (3.6) 

7T 

where 
and 

iNLO 



x / dyC^(2y- l)<f>(y,d), (3-7) 



(°) I a 

oo / / 2\ \ 7n /Pa x 

(x,/i|)=x(l-x)E / Prfr / dyC^(2y-l)<p(y^l) 



c3/ 2 (2x _ n (2^ + 3) / _ om\ hp & {0 ; 



oo 



+ J^J C " 2[2X ~ l) (k + l)(k + 2) SUA) C ™ 
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(3.8) 



with abbreviations 



' as(jJ> 



2^\l+(7£ 0) -7i 0) )//3o' 



(3.9) 



cff = (2n + 3) 



^ -Po + AC F A kn , 2C F (A fcn -^(A; + 2)+^(l)) 



(A;-n)(A; + n + 3) 



+ 



(n + l)(n + 2) 



A fcn = V( fe + 3 +4 ) - V>(^) + 2^(fc - n) - + 2) - ^(1) 



(3.10) 
(3.11) 



and X)' denoting the sum running only over even n, while C'^ 2 (z) are Gegenbauer poly- 
nomials of order 3/2 and Cf = 4/3. In Eqs. ( |3~7| - |3.10D , 7^ are the usual anomalous 
dimensions: 



7i 0) = C* 



3 + 



(ra + l)(ra + 2) 



n+1 1 ' 

-4E 



i=l 



(3.12) 



while (3q and /?i are the first two terms in the expansion of the QCD /^-function, with (3 
given by (|2.3j), and 



(3.13) 



fa. = 102- Y n f- 
The function if)(z) appearing in ( |3.10| ) and (|3.11|) is defined as 



ij(z) = -(\nT(z)). 



(3.14) 



For nj = 3 and for the values of n we are interested in, we have taken the values of the 
anomalous dimension from []25 , |26| : 



7? ) = . 



7^ = 111.03. 



7l 1} = 150.28 . 



(3.15) 



To proceed, we expand the distribution amplitude (j)(x, /jLq) in terms of the Gegenbauer 
polynomials C^ 2 (2x — 1) (the eigenfunctions of the LO evolution equation) 



(x) = 0(x, /4) = 6x(l - ar) £ 'B n C 3 J 2 (2x - 1) . 



(3.16) 



n=0 



Owing to the orthogonality of the Gegenbauer polynomials, the expansion coefficients 
B n are given by 



Br, 



4(2n + 3) f 1 



6 (n + l)(n + 2) 



/ dxCl ,2 (2x- l)Mx) 
Jo 



(3.17) 



and the normalization condition ( |3.2j ) then implies B = 1. Substituting ( |3.16j ) into 
gives 



lLO, 



x,&) = 6x(l - x) £ 'SJ°(a£) 2 (2* - 1) 



(3.18) 



n=0 
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where 



B n 



' ois(^o) 



(3.19) 



so that, obviously, B^° = 1. Next, by taking ( |3.16| ) into account (|3[ 



becomes 



where 



and 



lNLO 



x, = 6x(l - x) E (2x 



fc=2 



fc-2 



I), 



n=0 



(3.20) 



(3.21) 





(2fc + 3) (n + 2)(n + l) (1) 



(fc + l)(fc + 2) 2(2n + 3) 



(3.22) 

(3.23) 



For the purpose of our calculation, we use the following four candidate distribution 
amplitudes (shown in Fig. [5]) as nonperturbative inputs at the reference momentum scale 
Hi = (0.5 GeV) 2 : 

as (x) = (fras(x,nl) = 6sc(l - x) , (3.24a) 
<p cz {x) = <f> cz {x, fi 2 ) = 6x(l - x) [5(2x - l) 2 ] , (3.24b) 

P2 (x) = P2 (x, Ato) = M 1 - [-0.1821 + 5.91(2x - l) 2 ] , (3.24c) 

<j) P3 (x) = <pp 3 (x, i4) = 6x(l - x) [0.6016 - 4.659(2s - l) 2 + 15.52(2a; - l) 4 ] . (3.24d) 

Here as (x) is the asymptotic distribution amplitude and represents the solution of the 
evolution equation ( |3.4j ) for p? F — > 00. The double-hump-shaped distribution amplitudes 
fyczix) and 4>p 2 (x) and the three- hump-shaped distribution amplitude 4>p3(x) have been 
obtained using the method of QCD sum rules [^,^]. As Fig. |^ shows, these distribution 
amplitudes, unlike (f> as (x), are strongly end-point concentrated. In the limit fi F — ► 00, 
they reduce to the asymptotic form <fi as (x). 

The pion candidate distribution amplitudes given by ( |3.24| ) are of the general form 



<j>{x) = <P(x^l) = <p as (x) [1 + B 2 Cl /2 (2x - 1) + B 4 Cl /2 (2x - 1)] 
with the corresponding coefficients 



(3.25) 



1>as 


[x): 


B n 


= 


n > 2 






CZ 


[x): 


B n 


= 


n > 2 


B 2 


= 2/3 


!>P2 


[x): 


B n 


= 


n > 2 


B 2 


= 0.788 


!>P3 


[x): 


B n 


= 


n > 4 


B 2 


= 0.7582 



B A = 0.3941 



(3.26a) 
(3.26b) 
(3.26c) 
(3.26d) 
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FIG. 5. The four candidate pion distribution amplitudes denned by Eqs. ( 3.24| ), chosen as 
nonperturbative inputs at the reference scale //q = (0.5GeV) 2 . 



with Bq = 1. Now, according to ( 3.18 ) and ( |3.20 ), the LO and NLO parts of the 
distribution amplitude ( p\25| ) read 



^(x, /4) = <p as (x) 1 + B^iti?) Cr(?x - 1) + B^(j4) Cr(2x - 1 



■>LO/ , .2 



3/2, 



>LO/..2 



-,3/2, 



i NLO 



(x,n 2 F ) = <P as (x) [sf LO ( M |) Cl ,2 {2x - 1) + i?r o G4) Cl'\2x - 1) 



where 



+ £'< LO (/4)Q 3/2 (2*-i) 

fc=6 



Bo 



-50/81 



B A 



' <*s(t4) 



-364/405 



(3.27) 



(3.28) 



(3.29) 



and 



5 2 LO (/4m(/4) + Q2o(/4), 

5f°(^) P 4 (/4) + Q 4 oG4) + B^°(fi 2 F ) Q 42 (/4) , 
Q fe o(/4) + £ 2 L °(/4) Qfc 2 (/i|) + 5 4 LO (/i|) Q M (/4) , > 6 . (3.30) 

Note that <f) NLO (x, fi 2 F ) represents the infinite sum of Gegenbauer polynomials even 
though the (f>(x) distribution is described by a finite number of terms. Since Bj? de- 
creases with k, for the purpose of numerical calculation one can approximate 4> N (x, jMp) 
by neglecting higher-order Gegenbauer polynomials (k > 100). 

A summary of our results for the four candidate distribution amplitudes with the NLO 
evolutional corrections included is shown in Fig. The dash-dotted curves correspond 
to the distribution amplitudes at the reference point fi F = /i 2 , = (0.5 GeV) 2 . The dashed 
and solid curves represent the distribution amplitudes evoluted to fi F = (2 GeV) 2 , with 
the difference that the former includes only LO evolutional corrections, whereas the 



< LO (/4) = 
sf i0 (/4) = 
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FIG. 6. Evolution of the four candidate pion distribution amplitudes: <fi as (x, fx F ), 
4>cz(x, /x F ), 4>p2{x, n F ), and (f)ps(x, ftp), with the running coupling constant as{fx F ) and three 
active flavors. The dash-dotted curves correspond to the distribution amplitudes given by 



( 3.24 ) and taken as nonperturbative input at the reference momentum scale fx^ = (0.5 GeV) z 



The dashed curves correspond to the distribution amplitudes at the momentum scale fx F = (2 



GeV) 2 with the LO evolutional corrections included according to ( 3.27|) . The solid curves 
correspond to the distribution amplitudes at the momentum scale fxp = (2 GeV) 2 with the 
NLO evolutional corrections taken into account according to Eqs. (|3.6| ) and ( |3.27H3?2g| ). 



latter includes NLO evolutional corrections. As it is obvious from ( |3.26a| ) and ( |3.27| - 
3.28 ), although <p as shows no LO evolution, there are tiny effects of the NLO evolution. 
For other distributions, the LO evolution is significant, and even the NLO evolution is 
nonnegligible. Any distribution amplitude evolves asymptotically (i.e., for ix 2 F — ► oo) 
into <p as (x), but the higher \j? F is, the "slower" this approach becomes. 



IV. CHOOSING THE FACTORIZATION AND THE RENORMALIZATION 

SCALES 

In this section we discuss various possibilities of choosing the renormalization scale 
ix r and the factorization scale ixp appropriate for the process under consideration. 

Before starting the discussion, let us now, with the help of Eqs. (|1.1|), (|2.1| ), and (^.6|), 
write down the complete leading-twist NLO QCD expression for the pion electromagnetic 
form factor with the /xr and ixp dependence of all the terms explicitly indicated. 

Generally, for the NLO form factor we can write 
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FAQ 2 , /4 &) = Fl 0) (Q\ /4 /4) + F^iQ 2 , /4 A) ■ (4-1) 



The first term in (4.1) is the LO contribution and is given by 

F(°)(Q 2 ,/i^^)= fdx fdy a s (» R )$ LO *(y,» F )TP(x,y,Q 2 )$ LO (x,fi). (4.2) 
Jo Jo 

The second term in (|4.1| ) is the NLO contribution and can be written as 

F^ (Q 2 , /4 fi) = F^ (Q 2 , /4 /4) + Fl lb) (Q 2 , /4 A) , (4-3) 

where 

F™ = fdx fdy^^L 

JO JO 7T 

x $ LO * (y, &) if (x, y, Q 2 ) if (x, y, ^/Q 2 , ,4/Q 2 ) $ LO (x, , (4.4) 

is the contribution coming from the NLO correction to the hard-scattering amplitude, 
whereas 

1 '"' a s{fi R )ocs{fi 



Fi lb \Q 2 ,,4/4) = f dx { dy 
Jo Jo 



F) 



x [s^fo a*|.) Tg\x, y, Q 2 ) $ LO (x, &) + /£) Tf(x, y, Q 2 ) $™(x, 

(4.5) 

is the contribution arising from the inclusion of the NLO evolution of the distribution 
amplitude. Now, if ( |2.4| ) is taken into account, the expression for the LO contribution 
of Eq. (|4.2j) can be written in the form 



F r(Q , hr, hf) = g * ^ y o rfx y o «s(mh) — - - — • ( 4 -6) 

Next, before going on to perform the remaining x and y integrations, we have to 
choose the renormalization scale fiR and the factorization scale y,p. In doing this, how- 
ever, there is considerable freedom involved. 

If calculated to all orders in perturbation theory, the physical pion form factor F n (Q 2 ), 
represented at the sufficiently high Q 2 by the factorization formula (|1.1|), would be in- 
dependent of the renormalization and factorization scales, fi R and /if, and so they are 
arbitrary parameters. Truncation of the perturbative expansion of F^iQ 2 ) at any finite 
order causes a residual dependence on these scales. Although the best choice for these 
scales remains an open question (the scale ambiguity problem), one would like to choose 
them in such a way that they are of order of some physical scale in the problem, and, 
at the same time, to reduce the size of higher-order corrections as much as possible. 
Choosing any specific value for these scales leads to a theoretical uncertainty of the 
perturbative result. 

In our calculation we approximate F n (Q 2 ) only by two terms of the perturbative 
series and hope that we can minimize higher-order corrections by a suitable choice of 
Hr and jjiF-, so that the LO term F^ (Q 2 , (j? R , fj? F ) gives a good approximation to the 
complete sum F^iiQ 2 ). 
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It should be noted that in F^(Q 2 , fi R , fj, F ), given by ( |4.2|) , fi R appears only through 
as((J? R ), whereas [if enters into the distribution amplitude (p(x, /j> f ). In F^(Q 2 , fj, R , jjp), 
given by fl4.3| - f4~5l) , a logarithmic dependence on the scales fi R and /if appears also 
through Tjj\x, y, fi 2 R /Q 2 , fi 2 F /Q 2 ). As seen from (|2.21|) and (|2.20|) , this dependence is 
contained in the terms 



fuv(x,y,Q 2 /(A 



fm(x, y, Q 2 / H 2 F ) = - [3 + \n(x y )] 



--ln(xy) + lngf 

i - . / . in, 



2 \n(xy)-\n ^ 2 



(4.7) 
(4.8) 



Being independent of each other, the scales fi R and fiF can be expressed in terms of Q 2 
as 



/4 

/4 



a(x,y) Q 2 

K x iv)Q 2 : 



(4.9) 
(4.10) 



where a(x,y) and b(x,y) are some linear functions of the dimensionless variables x and 
y (quark longitudinal momentum fractions). 

We next discuss various possibilities of choosing the scales fiR and \ip separately. The 
simplest and widely used choice for the scale \i R is 



/4 



Q 2 



(4.11) 



the justification for the use of which is mainly pragmatic. 

Physically, however, the more appropriate choice for fijz, would be that one correspond- 
ing to the characteristic virtualities of the particles in the parton subprocess, which is 
considerably lower than the overall momentum transfer Q 2 (i.e., virtuality of the probing 
photon). 

It follows from Figs. [I] and |2| that the virtualities of the gluon line (line 4) and 
the internal quark line (line 2) of diagram A in Fig. |2| are given by xyQ 2 and yQ 2 , 
respectively. Now, if instead of using ( [4.11| ) we choose to be equal to the gluon 
virtuality, i.e., 



/4 



xyQ 2 



(4.12) 



then the logarithmic terms in (|4.7j ) vanish. 

As it is well known, unlike in an Abelian theory (e.g., QED), where the effective 
coupling is entirely renormalized by the corrections of the vector particle propagator, in 
QCD the coupling is renormalized not only by the gluon propagator, but also by the 
quark-gluon vertex and quark-propagator corrections. It is thus possible to choose n 2 R 
as the geometrical mean of the gluon and quark virtualities ||: 



/4= ^(xyQ 2 )(yQ 2 ) 
Alternatively, we can make a choice 

V 2 R = e~^xyQ 2 , 



(4.13) 



(4.14) 
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as a result of which the function fuv, 



T^ J (x,y, fJl 2 R /Q 2 ,fJ 2 F /Q 2 ) defined by ( ^20|) becomes n 



given by ( |4.7| ), vanishes identically. In this case, 
f independent. This is an example 
of choosing the renormalization scale according to the Brodsky-Lepage-Mackenzie (BLM) 
procedure In this procedure, the renormalization scale fi R best suited to a particular 
process in a given order can be determined by computing vacuum-polarization insertions 
in the diagrams of that order. The essence of the BLM procedure is that all vacuum- 
polarization effects from the QCD (3 function are resummed into the running coupling 
constant. 

Let us just mention at this point that in addition to the BLM procedure, two more 
renormalization scale-setting procedures in PQCD have been proposed in the literature: 
the principle of fastest apparent convergence (FAC) [^9|, and the principle of minimal 
sensitivity (PMS) [30] . The application of those three quite distinct methods can give 



strikingly different results in practical calculations |3T . 

As for the factorization scale [i F , it basically determines how much of the collinear 
term given in ( |4.8| ), is absorbed into the distribution amplitude. A natural choice for 
this scale would be 



/4 = Q 2 , (4.15) 
which eliminates the logarithms of Q 2 /fi 2 F . More preferable to ( }4.15|) is the choice 

/i| = vWg 2 , (4.16) 
, vanish. 



where the coupling constants «s(/^) and as(n F ) appear 



which makes the function fj R , given by 

A glance at Eqs. (|4lH4 
under the integral sign, reveals that any of the choices of fi R given by ( f4.12| - [L"T4| ), and 
the choice of \i F given by Q4.16| ), leads immediately to the problem if the usual one- loop 
formula fl2.2p for the effective QCD running coupling constant is employed. Namely, 
then, regardless of how large Q 2 is, the integration of Eqs. (|4~l"H4.5|) allows a^/z^) to be 
evaluated near zero momentum transfer. Two approaches are possible to circumvent this 
problem. First, one can choose \j? r and fi 2 F to be effective constants by taking fi 2 R = (fi R ) 
and n 2 F = (fi 2 F ), respectively. Second, one can introduce a cutoff in formula Q2.2| ) with 
the aim of preventing the effective coupling from becoming infinite for vanishing gluon 
momenta. 

If the first approach is taken, Eqs. (^4.12| - f3~T4"D and ( |4.16|) get replaced by the averages 



/4 



(xyQ' 



and 



(xyQ 2 ) (yQ 2 ) 



-5/3. 



(J>F 



xyQ' 



(xy)Q 2 



respectively. Taking into account the fact that (x y ) = (x ) (y } and (x ) 
possible to write Eqs. ( ^.17| ) and ( (4.18[ ) in the respective forms: 



(4.17a) 

(4.17b) 
(4.17c) 

(4.18) 
(y), it is 
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fx 2 R =(x) 2 Q 2 , (4.19a) 

V 2 R =(xf 2 Q 2 , (4.19b) 

fi 2 R = e- 5 / 3 (x) 2 Q 2 , (4.19c) 

and 

li 2 F = (x)Q 2 . (4.20) 

The key quantity in the above considerations is (x), the average value of the mo- 
mentum fraction. It depends on the form of the distribution amplitude, and there is no 
unique way of defining it. A possible definition is 

/ dx(f)(x, n%)x 

<z)G4) = JjL n • ( 4 - 21 ) 

dx 4>(x, fx 2 F ) 



o 



Owing to the fact that all distribution amplitudes under consideration are centered 
around the value x = 0.5, it follows trivially from (|4.21|) that for all of them 

(x)(/4) = 0.5. (4.22) 

An alternative way of defining (x), motivated by the form of the LO expression for the 
pion form factor ( |4.6[ ) , is 

/•i LO (x,/4) _ 

/ dx — x i 

{w } = V ix Lo ^,^) - = w+wm+wm+^} ■ (4,23) 

Jo X 

It should be noted, however, that this formula can be straightforwardly applied only 
if fi 2 F = Q 2 . On the other hand, if instead of fip = Q 2 one chooses the factorization 
scale to be as given by Q4.20|) , then Eqs. (|4.20|) and (|4.23|) form a nontrivial system 



of simultaneous equations. According to ( |4.23| ), one obtains (x) as (Q 2 ) = 1/3 for any 
Q 2 , while 0.242 < (x) cz (Q 2 ) < 0.262 for 4GeV 2 < Q 2 < lOOGeV 2 (similar values are 
obtained for <pp 2 and (pps). 

When using the <fi as (x, (ip) distribution it appears reasonable to take (x) (fi 2 F ) = 
(x) as =l/2. This can be justified on the grounds that this distribution is concentrated 
around x = 0.5, and is characterized by a very weak evolution. On the other hand, for 
the end-point concentrated distributions 4>cz(x, fip), 4>p2{x,n 2 F ) and <fip3(x, /i 2 ^), which 
exibit sizable evolutional effects, it is more appropriate to take (x)(n 2 F ) as given by Eq. 

(US). 



As stated above, the divergence of the effective QCD coupling as(/z|,), as given by 
72]), is the reason that it is not possible to use the choices of fi 2 R given by (|4. 12| - fi~T4l) 



and fj, F given by ( f4.16p . Equation fl2T2] ) does not represent the nonperturbative behavior 
of as^p) for small /zfj, and a number of proposals have been suggested for the form of 
the coupling constant in this regime |B2|. The most exploited parameterization of the 
effective QCD coupling constant at low energies has the form 
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a s{» 2 R ) 



ft In 



Att 

A 2 



(4.24) 



where the constant C encodes the nonperturbative dynamics and is usually interpreted 
as an "effective dynamical gluon mass" m g . For fi R ^> C 2 , the effective coupling in 
d4.24j ) coincides with the one-loop formula (|2.2| ), whereas at low momentum transfer this 
formula "freezes" to some finite but not necessarily small value. 



In view of the confinement phenomenon, the modification (|4.24j) is very natural: the 
lower bound on the particle momenta is set by the inverse of the confinement radius. 
This is equivalent to a strong suppression for the propagation of particles with small 
momenta. Thus, in a consistent calculation in which ( f4.24| ) is assumed a modified gluon 
propagator should be used: 



1 



k 2 



(4.25) 



However, if one attempts to calculate the LO prediction for the pion form factor making 
use of (|4.24j) and fl4.25|) , and taking m g = 300 — 500 MeV, the result obtained is by a 
factor of 10 lower than the experimental value, questioning the applicability of such an 
approach [|l^ . 

It has recently been shown that spontaneous chiral symmetry breaking imposes rather 



a severe restriction on the idea of freezing ||33f| . The authors of |33[ argue that before any 
argument based on a particular form of the freezing coupling constant is put forward, one 
should check that the dynamical origin (mechanism) of the freezing is such that enough 
chiral symmetry breaking can be produced. 

Considering the discussion above, calculations with the frozen coupling constant seem 
to need a more refined treatment and will not be considered in this paper. 



V. COMPLETE NLO NUMERICAL PREDICTIONS FOR THE PION FORM 

FACTOR 



Having obtained all the necessary ingredients in the preceding sections, now we put 
them together and obtain complete leading-twist NLO QCD numerical predictions for 
the pion form factor. For a fixed distribution amplitude, we analyze the dependence of 
our results on the choice of the scales fi R and /ip. 

By inserting ( |3.18|) and ( |3.20|) into ( [4. 2D and Q1.5D , taking into account (|2~4l ), taking 
the scales fi 2 R and up to be effective constants (as explained in Sec. IV), and performing 
the x and y integration, we find that Eqs. fl4.2|) and ( [4.5|) take the form 



and 



\n=0 



Q 2 ^ lb) (Q 2 ,/4/4) = lQflas^ 2 R )as^ 2 F ) £ 'B 



\n=0 




k=2 



(5.1) 



(5.2) 



21 



For a distribution amplitude of the form given in (|3.25| - |3.30|) , the above expression 
reduces to 



Q 2 4 lb) (Q 2 ,»%» 2 F ) = lQfla s ^ R )a s ^ 2 F ) \l + B^^ 2 F ) + B L A °^ 2 F ) 



k=6 



(5.3) 



(5.4) 



As for the part of the NLO contribution arising from the NLO correction to the hard- 
scattering amplitude, by inserting ( |2.20| ) and ( |3.27| ) into ( f4.4| ) and performing the x and 
y integration one obtains 



~\l + B™(n F ) + B™(^ 2 



+- 



25 



B L 2 °{» 



91 
15 



l + B^{n F ) + B^^ F 



LOi 



111 



/4 
Q 2 



+6.58 + 24.99 B^°(n 2 F ) + 21.43 {B^°(^ F )) 2 

+32.81 B%°(fi F ) + 32.55 (5 4 LO (/4)) 2 + 53.37 5 2 LO (/4) 5 4 LO (/4)} 

Then the NLO contribution to the "almost scaling" combination Q 2 F n (Q 2 ) is given by 

Q 2 F«(Q 2 ,^,/,|) = Q 2 F^\Q 2 ^ 2 R ^ 2 F ) + Q 2 F^\Q 2 ^%&) , (5.6) 



while the total NLO prediction reads 

Q 2 F n {Q 2 ,fi 2 R ,fi 2 F ) = Q 2 F?\Q 2 ^ 2 R ^ 2 F ) + Q 2 F£\Q 2 ^ 2 R ^ 2 F ) . 



(5.7) 



For the purpose of this calculation we adopt the criteria according to which a per- 
turbative prediction for F n (Q 2 ) is considered reliable provided the following two require- 
ments are met: first, corrections to the LO prediction are reasonably small (< 30%); 
second, the expansion parameter (effective QCD coupling constant) is acceptably small 
(as(fi 2 R ) < 0.3 or 0.5). Of course, one more requirement should be added to the above 
ones: consistency with experimental data. This requirement, however, is not of much 
use here since reliable experimental data for the pion form factor exist for Q 2 < 4 GeV 2 , 
i.e., well outside the region in which the perturbative treatment based on Eq. ( |1.1| ) is 
justified. 

Currently available experimental data for the spacelike pion electromagnetic form 
factor F 1X {Q 2 ) are shown in Fig. [F[ The data are taken from Bebek et al. |35| and 
Amendolia et al. [BB]. As stated in |35[|, the measurements corresponding to Q 2 = 



6.30GeV 2 and Q 2 = 9.77GeV 2 are somewhat questionable. Thus, effectively, the data 
for F n (Q 2 ) exist only for Q 2 in the range Q 2 < 4GeV 2 . The results of Ref. ||3l| were 
obtained from the extrapolation of the 7*p — > TT + n electroproduction data to the pion 



pole. It should also be mentioned that the analysis of Ref. |35| is subject to criticism 
which questions whether F n (Q 2 ) was truly determined for 1 < Q 2 < 4GeV 2 |37j (but 
see also Ref. [[38|). The new data in this energy region are expected from the CEBAF 
experiment E-93-021. 
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FIG. 7. Presently available experimental data for the spacelike pion electromagnetic form 
factor. 



A. Predictions obtained with p? R 



/4 



Q 2 



The first NLO prediction for the pion form factor was obtained in Ref. 
the MS renormalization scheme and the choice fi R 



f4 



. Using 

Q 2 it was found that for 



distribution (with the evolution of the distribution amplitude neglected), the 



the <t) as (x] 
perturbative series took the form 

Q 2 F n (Q 2 ) = (0.43 GeV 2 ) a s (Q 2 ) [1 + 2.10 a s (Q 2 



(5.1 



which is in agreement with our result given by (|5.7|) , (|5.3|) , and ( |5.5|) . The conclusion 
based on this result was that a reliable result for F n (Q 2 ) was not obtained until as{Q 2 ) ~ 
0.1, or with Ajj-g = 0.5 GeV, Q 2 = 10000 GeV 2 . This prediction has been widely cited 
in the literature and initiated a lot of discussion regarding the applicability of PQCD to 
the calculation of exclusive processes at large momentum transfer. With the presently 
accepted value of 0.2 GeV for Aj^g, we find that the criteria from Ref. || are satisfied 
for Q 2 ~ 1600 GeV 2 , and that for Q 2 = 700 GeV 2 , the corrections to the LO prediction 
are of order 30%. Thus, this result shows that for the choice of the renormalization and 
factorization scales fi R = fi 2 F = Q 2 , the region in which perturbative predictions can be 
considered reliable is still well beyond the region in which experimental data exist. The 
inclusion of the distribution amplitude evolution effects, although extremely important 
for the end-point concentrated amplitudes, does not change this conclusion. 

Numerical results of our complete NLO QCD calculation, obtained using the four 



candidate distribution amplitudes, with fj, R 



/4 



Q 2 , and for Q 2 > 4 GeV 2 



arc 



displayed in Tables [T| |T|, [TV], and |V|. The entries in these tables include various contri- 
butions given by Eqs. Q5.3| - |5.71 ), comprising the full NLO result. A comparison of our 
results with presently available experimental data is shown in Fig. |8|. The ratio of the 
NLO to the LO contribution to i^(Q 2 ), i.e., Fjp (Q 2 ) / (Q 2 ) , as a useful measure of 
the importance of the NLO corrections, is plotted as a function of Q 2 in Fig. []. The 
shaded areas appearing in Figs. || and |9| denote the region of Q 2 where as(Q 2 ) > 0.3. 
We take that outside of this region the effective coupling is acceptably small. 
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TABLE II. Complete leading-twist NLO QCD results for the pion form factor, Q 2 F 7T (Q 2 ), 
obtained using the (pas(x, fJ-'p) distribution amplitude and assuming /i 2 ^ = p? F = Q 2 . 



Q 2 a 3 (p 2 R ) Q 2 Fl '(Q 2 ) Q 2 Fl '(Q 2 ) Q 2 Fr\Q 2 ) Q 2 F^(Q 2 ) F?\q 2 ) / F^ (Q 2 ) Q 2 F 7V (Q 2 ) 



[GeV 2 ] 




[GeV 2 ] 


[GeV 2 ] 


[GeV 2 ] 


[GeV 2 ] 


% 


[GeV 2 


4 


0.303 


0.131 


0.083 


-0.001 


0.082 


62.8 


0.213 


6 


0.279 


0.120 


0.070 


-0.001 


0.069 


57.6 


0.189 


8 


0.264 


0.114 


0.063 


-0.001 


0.062 


54.4 


0.176 


10 


0.253 


0.109 


0.058 


-0.001 


0.057 


52.2 


0.166 


20 


0.225 


0.097 


0.046 


-0.001 


0.045 


46.2 


0.142 


30 


0.211 


0.091 


0.040 


-0.001 


0.039 


43.3 


0.130 


40 


0.202 


0.087 


0.037 


-0.001 


0.036 


41.4 


0.123 


50 


0.196 


0.085 


0.035 


-0.001 


0.034 


40.1 


0.118 


75 


0.185 


0.080 


0.031 


-0.001 


0.030 


37.9 


0.110 


100 


0.178 


0.077 


0.029 


-0.001 


0.028 


36.4 


0.105 



TABLE III. Same as Table || but for the 4>cz(x, [i 2 F ) distribution amplitude. 



Q 2 


"SOD 


Q 2 Fi 0) (Q 2 ) 


Q 2 Fi la) (Q 2 ) 


Q 2 Fl U) (Q 2 ) 


Q 2 F^(Q 2 ) 


Fi 1] (Q 2 )/Fi°\Q 2 ) 


Q 2 F 7T (Q 


[GeV 2 ] 




[GeV 2 ] 


[GeV 2 ] 


[GeV 2 ] 


[GeV 2 ] 


% 


[GeV 2 ] 


4 


0.303 


0.248 


0.241 


-0.015 


0.225 


90.8 


0.474 


6 


0.279 


0.222 


0.195 


-0.014 


0.181 


81.6 


0.403 


8 


0.264 


0.206 


0.170 


-0.013 


0.157 


76.1 


0.363 


10 


0.253 


0.195 


0.153 


-0.012 


0.141 


72.2 


0.336 


20 


0.225 


0.167 


0.115 


-0.011 


0.104 


62.2 


0.271 


30 


0.211 


0.154 


0.098 


-0.010 


0.088 


57.4 


0.243 


40 


0.202 


0.146 


0.089 


-0.009 


0.079 


54.3 


0.225 


50 


0.196 


0.140 


0.082 


-0.009 


0.073 


52.2 


0.213 


75 


0.185 


0.131 


0.072 


-0.008 


0.064 


48.7 


0.194 


100 


0.178 


0.125 


0.065 


-0.008 


0.058 


46.4 


0.182 



TABLE IV. Same as Table but for the 4>p2{x, fi F ) distribution amplitude. 



Q 2 




) Q 2 F^{Q 2 ) Q 2 Fi la) (Q 2 ) 


Q 2 Fi lb) (Q 2 ) 


q 2 fP(q 2 ) 


F^(Q 2 )/Fi°\Q 2 ) 


Q 2 F 7T (Q 


[GeV 2 ] 




[GeV 2 ] 


[GeV 2 ] 


[GeV 2 ] 


[GeV 2 ] 


% 


[GeV 2 ] 


4 


0.303 


0.274 


0.278 


-0.019 


0.259 


94.6 


0.532 


6 


0.279 


0.244 


0.224 


-0.017 


0.207 


85.0 


0.451 


8 


0.264 


0.226 


0.195 


-0.016 


0.179 


79.1 


0.404 


10 


0.253 


0.213 


0.175 


-0.015 


0.160 


75.0 


0.374 


20 


0.225 


0.182 


0.130 


-0.013 


0.117 


64.5 


0.300 


30 


0.211 


0.167 


0.111 


-0.012 


0.100 


59.4 


0.267 


40 


0.202 


0.158 


0.100 


-0.011 


0.089 


56.2 


0.247 


50 


0.196 


0.152 


0.093 


-0.011 


0.082 


54.0 


0.234 


75 


0.185 


0.141 


0.081 


-0.010 


0.071 


50.2 


0.212 


100 


0.178 


0.134 


0.074 


-0.009 


0.064 


47.9 


0.199 
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FIG. 8. Comparison of the complete leading-twist NLO QCD predictions for the pion form 
factor, Q 2 Ft T (Q 2 ), obtained using the four candidate distribution amplitudes, with the presently 
available experimental data. The shaded area denotes the region of Q 2 in which as(Q 2 ) >0.3. 
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FIG. 9. The ratio of the NLO to the LO contributions to the pion form factor, 
Fn (Q 2 ) / Fn (Q 2 ) , for the four candidate distribution amplitudes. The shaded area denotes 
the region of Q 2 in which as(Q 2 ) >0.3. 
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TABLE V. Same as Table || but for the (j)pz{x, /j, 2 ?) distribution amplitude. 



Q 2 


" Vr"it 


) q 2 fP(q 2 ) 


Q 2 Fi la) (Q 2 ) 


Q 2 Fi lb) (Q 2 ) 


Q 2 Fi 1) (Q 2 ) 


fP(q 2 )/f^ 

x TT \ y °C J j ■*• TV 


(Q 2 ) Q 2 FJQ 


[GeV 2 l 




[GeV 2 l 


[GeV 2 l 


fGeV 2 l 


[GeV 2 l 


% 


fGeV 2 l 


A 

-± 


u.ouo 


u.ooo 


A09 


09Q 


379 


1110 


70S 
U. 1 uo 


a 
u 


D 97Q 


9Q^ 


318 


n D9fi 

-u.uzu 


9Q9 
u.zyz 


1 00 


^R7 

U.OO I 


Q 





979 


D 973 




9zLQ 


Q1 S 


^91 

U.OZ 1 


1 




n 9^ 


9zLA 


n 093 


999 


ou. / 


A77 






n 9i ^ 


n 1 7s 

U. 1 / o 


n ni q 


1 


73 S 


373 
u.o / o 


30 


n 01 1 

0.211 


0.19b 


0.150 


-0.017 


0.133 


67.7 


0.328 


40 


0.202 


0.184 


0.134 


-0.016 


0.118 


63.8 


0.302 


50 


0.196 


0.176 


0.123 


-0.015 


0.108 


61.1 


0.284 


75 


0.185 


0.163 


0.106 


-0.014 


0.092 


56.6 


0.255 


100 


0.178 


0.154 


0.096 


-0.013 


0.083 


53.7 


0.237 



It is evident from Figs. [5] and [5] and Tables fJ^V] that the leading-twist NLO results 
for the pion form factor obtained with fi 2 R = fi 2 F = Q 2 display the following general 
features. First, the results are quite sensitive to the assumed form of the pion distribution 
amplitude. Thus, the more end-point concentrated distribution amplitude is, the larger 
result for the pion form factor is obtained, and also the NLO corrections are larger 
(which has already been obvious looking at Eqs. ( |5.3| - |5.7| )). Second, whereas the NLO 
correction arising from the corrections to the hard-scattering amplitude are positive, 
the corrections due to the inclusion of the evolutional corrections to the distribution 
amplitude are negative, with the former being generally an order of magnitude larger 
than the latter. Thus, in all the cases considered, the full NLO correction to the pion 
form factor is positive, i.e., its inclusion increases the LO prediction. 

We now briefly comment on the results obtained with each of the four distribution 
amplitudes. 

Table |TJ, which corresponds to the <fi as (x,Q 2 ) distribution amplitude, shows that 
the NLO correction Q 2 F«(Q 2 ) is rather large (36% at Q 2 = 100 GeV 2 ). Most of the 
contribution to Q 2 Fj T 1 \Q 2 ) is due to the NLO correction to the hard-scattering ampli- 
tude Q 2 Fj la )(Q 2 ), while the contribution Q 2 F^ lb \Q 2 ) arising; from the NLO evolutional 
correction of the distribution amplitude is rather small, being of order 1%. The ratio 
F^\Q 2 ) I F^\Q 2 ) « 30% is reached at Q 2 » 500 GeV 2 . 

The results derived from the (pcz(x, Q 2 ) distribution are presented in Table |TTT|. The 
full NLO correction Q 2 F^(Q 2 ) is somewhat larger than for (p a s(x,Q 2 ) (and at Q 2 = 
100 GeV 2 it amounts to 46%). The ratio (Q 2 ) / (Q 2 ) is greater than 30% until 
Q 2 « 2400 GeV 2 . It is important to observe that the evolutional corrections, especially 
the LO ones, are rather significant in this case. Also, as it is seen from Table |TT], the 
NLO evolutional correction Q 2 F^ lb \Q 2 ) is of order »s 6%, i.e., nonnegligible. To show 
the importance of the correction arising from the inclusion of the distribution amplitude 
evolution, the results for Q 2 F n (Q 2 ) and the ratio F^ (Q 2 ) / F^ (Q 2 ) obtained using the 
4>cz(%), (fic°(x,Q 2 ), and (ficz(%, Q 2 ) distributions, are exhibited in Figs. |TU| and [TT], 
respectively. 

The results based on the <fip2(x,Q 2 ) and 4>P3(x,fip) distributions are listed in Table 
|TV| and |V|, respectively. As it can be easily seen by looking at Figs. [8] and |9| and by 
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FIG. 10. Leading-twist NLO QCD predictions for the pion form factor, Q 2 F 7T (Q 2 ), obtained 
with the (f>cz{x), 4>cz( x i Q 2 )i an d <frcz(x, Q 2 ) distributions. The shaded area denotes the region 
of Q 2 in which a s (Q 2 )> 0.3. 
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FIG. 11. The ratio of the NLO to the LO contributions to the pion form factor, 
Fi l) (Q 2 )/Fi°\Q 2 ), obtained for the 4>cz(x), (j)^{x,Q 2 ), and (j) C z{x,Q 2 ) distributions. The 
shaded area denotes the region of Q 2 in which as(Q 2 ) >0.3. 
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comparing the corresponding entries in Tables |(| and [IV], the results obtained with 
the 4>p2(x, Q 2 ) and (f>cz(x, Q 2 ) distributions are practically the same qualitatively, while 
differ quantitatively by a few percent. From Table [V] and Figs. |8] and |, one can see 
that the behavior of the results obtained with the (f>P3(x, fi 2 F ) distribution is qualitatively 
similar to the behavior of the results obtained with other two end-point concentrated 
distributions. For this reason, we leave the 4>p2{x,Q 2 ) and (f)p3(x, ii 2 F ) distributions out 
of our further consideration. 

In view of what has been said above, we may conclude the following. If the pion is 
modeled by the 4> as (x, fip), (f)cz(x, fi 2 F ), <fip2(x, /^f), or <pp3(x, n 2 F ) distribution amplitude, 
and if the renormalization and factorization scale are chosen to be pL R = fi 2 F = Q 2 , one 
finds that the NLO corrections to the lowest-order prediction for the pion form factor are 
large. The NLO predictions obtained cannot be made reliable, i.e., Fjp (Q 2 ) / F^ (Q 2 ) less 
than, say, 30 %, until the momentum transfer Q ^> 10 GeV is reached. Based on these 
findings and considering the region of Q 2 in which the data exist, it is clear that we are not 
in a position to rule out any of the four distributions considered. One can only note that 
the predictions for Q 2 F n (Q 2 , n R , ii 2 F ) obtained with the (p as (x, fj, F ) distribution are below 
the trend indicated by the existing experimental data, while the end-point concentrated 
distributions (ficz{%, Hf), 4 ) P2{x, Hp), and (pp 3 (x, fi 2 F ) give higher predictions. It is worth 
mentioning here that the theoretical predictions for the photon-to-pion transition form 
factor F^IQ 2 ) are in very good agreement with the data, assuming the pion distribution 
amplitude is close to the asymptotic one, i.e., (p as (x, fi 2 F ) [j39fl . 



B. Predictions obtained using ^jl 2 r = aQ 2 and fip = bQ 2 

In this subsection we present a detailed analysis of the dependence of the complete 
leading— twist NLO predictions for the pion form factor obtained with the <f) a s{x, Hf) and 
(ficzix, Hf) distributions, on the renormalization and factorization scales, /i# and fip- In 
the following we shall restrict our attention to two most exploited pion distribution 
amplitudes 4> a s(x,fj, F ) and 4>cz(x, fi F ). 

In the preceding subsection, we have found that the NLO corrections calculated using 
these two distributions are large, especially for the latter. The reason for this lies in the 
fact that the renormalization scale choice fi 2 R = Q 2 is not appropriate one. Namely, 
owing to the partitioning of the overall momentum transfer Q 2 among the particles in 
the parton subprocess, the essential virtualities of the particles are smaller than Q 2 , so 
that the physical renormalization scale, better suited for analyzing the process under 
consideration, is inevitably lower than the external scale Q 2 . 

Compared with the scale fip, the factorization scale \i F turns out to be of secondary 
importance. 

A characteristic feature of the asymptotic distribution as ( that in the LO it 

shows no evolution. A consequence of this, as it can be seen from Eqs. ( |5.3| - |5.7|) , is that 



the NLO predictions for the pion form factor based on this distribution are essentially 
independent of the factorization scale fiF- Namely, the only dependence on this scale is 
contained in the term arising from the NLO evolutional effects, which are tiny, as we have 
seen in Sec. III. On the other hand, the predictions calculated with the 4>cz{%, 1^%) are 
^f dependent, but this dependence turns out to be very weak. Thus, for a given value 
of fip, variation of the value of /ip in the range Q 2 / '4 < ftp < Q 2 leads to practically the 
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FIG. 12. The numerical results for Q 2 F^ (Q 2 ) obtained with the 4> a s(x, fip) distribution 
amplitude and the choices of fj, R given by Eqs. (14.19 ), with /xp = Q 2 and (x) as = 1/2. The 
solid curve (included for comparison) corresponds to the case fj, R = /J,p = Q 2 considered in the 
preceding subsection. The shaded area denotes the range of the LO prediction Q 2 F^ ] {Q 2 ) for 
/4/Q 2 € [e- 5 /3^) 2 ,l . 



same results. Therefore, when using the </> as (x, (i 2 F ) and (f>cz{x, Hf) distributions, one is 
allowed to set jjb 2 F = Q 2 , for all practical purposes. 



1. Examining the renormalization scale dependence of the NLO corrections 



The three specific physically motivated choices of (j,r, given by (|4.19| ), can be conve- 
niently written as 

ji 2 R = a Q 2 , (5.9a) 

where 

a G {(x) 3/2 ,(x) 2 ,e- 5/3 (x) 2 }. (5.9b) 

Our complete leading-twist NLO predictions for the pion form factor calculated with 
the 4>as{%, Hf) distribution and with the three specific values of the renormalization scale 
are summarized in Figs. [L3| and [14], showing the results for Q 2 F T (Q 2 ) and the ratio 



F^ (Q 2 ) I F( \Q 2 ) , respectively. For the purpose of the discussion of the effects of the 



inclusion of the NLO corrections on the LO predictions, we have also included Fig. |12 
showing the LO predictions for the three different values of fiR given in Eq. ( |5.9|) . The 
corresponding results but obtained assuming the 4>cz(x,fip) distribution are displayed 



in Figs. M, M, and 17 



The solid curves in Figs. [12] - [L7| correspond to the results with n 2 R = n 2 F = Q 2 
obtained in the preceding subsection, are included for comparison. The other curves 
refer to the choices of n 2 R given by (|5.9j). 
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FIG. 13. Leading-twist NLO QCD results for Q 2 F V {Q 2 ) obtained with the cf) as {x, ji 2 F ) dis- 
tribution amplitude and the choices of /j, r given by Eqs. Q4.1S| ), with p? F = Q 2 and (x) as = 1/2. 
The solid curve (included for comparison) corresponds to the = Q 2 considered 

in the preceding subsection. The shaded area denotes the range of the total NLO prediction 
Q 2 F 7T (Q 2 ) where the upper limit corresponds to fi 2 R = ^ 2 PMS obtained from (|5.11| ). 
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FIG. 14. The ratio F^(Q 2 ) / F^\q 2 ) obtained with the 



n%) distribution amplitude 



and the choices of fi 2 R given by Eqs. (|4.19|) , with fi 2 F = Q 2 and (x) as = 1/2. The solid 
curve (included for comparison) is for the result corresponding to n 2 R = up = Q 2 obtained 
in the preceding subsection, and the shaded area denotes the region of the predictions which 



correspond to a 5 (/x|)<0.5 (a s (/^) <0.3), while \F^ l) (Q 2 ) / FP (Q 2 )\ < 0.3 
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FIG. 15. The numerical results for Q 2 F^ (Q z ) obtained with the <j)cz(x, [ip) distribution 
amplitude and the choices of n R given by Eqs. (f4.19| ), with [i F = Q 2 and (x) cz (Q 2 ) calcu- 
lated according to ( [4.23 ). The solid curve (included for comparison) corresponds to the case 
Q 2 considered in the preceding subsection. The shaded area denotes the range of 



/4 



the LO prediction Q 2 Fi '(Q 2 ) for \i\jQ 2 G e" 5 / 3 (x) 2 , 1 



Some general comments concerning the results presented above are in order. First of 
all, it is interesting to note that the predictions based on the (f> as (x, Hf) an d (pcz(x, 
distributions amplitudes almost share the qualitative features, although these two dis- 
tributions are quite different in shape. Next, by looking at Figs. [13] and |l^, one notices 
that the full NLO result for Q 2 F n (Q 2 ), shows a very weak dependence on the value of 
fx R , and that it increases with decreasing fj, R . As far as the ratio F^> (Q 2 ) / F^ ' (Q 2 ) is 



concerned, as evident from Figs. [14] and [T7|, the situation is quite different. It is rather 



sensitive to the variation of fi R , and decreases as fi 2 R decreases. The choice fi 2 R = Q 2 , rep- 
resented by solid curves, when compared with the other possibilities considered, leads to 
the lowest value for Q 2 F n (Q 2 ), and to the highest value for the ratio F^ (Q 2 ) / F^ (Q 2 ) . 
In contrast to that, the choice of the BLM scale, represented by the dashed-dotted curve, 
leads to somewhat higher values of Q 2 F n (Q 2 ), but to considerably lower values of the 
ratio F«(Q 2 )/Fi°)(Q 2 ). 

Adopting the previously stated criteria, we now comment on the reliability of the 
NLO predictions for Q 2 F n (Q 2 ) displayed in Figs. [13] and [16]. Imposing the requirements 
F^ l \Q 2 )/F^(Q 2 )\ < 0.3, and a s {^ 2 R ) < 0.3, we find from Fig. [I| that for the os (x, pip) 

distribution the results corresponding to a = (x ) 3 ^ 2 , (x ) 2 , and e~ 5 ^ 3 (x} 2 , become reliable 
for the momentum transfer Q 2 > 50, 25, and 90 GeV 2 , respectively. It might be that the 
requirement as(fi 2 R ) < 0.3 is too stringent. Thus, by relaxing it and taking as(n R ) < 0.5 
instead, we find that the results corresponding to a = (x) 3//2 , (x) 2 , and e~ 5//3 (a;) 2 , 
become reliable for the momentum transfer Q 2 > 50, 25, and 15 GeV 2 , respectively. 

Applying the same criteria to the case of the (pcz(x, fi 2 R ) distribution, we find from 
Fig. [17] that the results for Q 2 F W (Q 2 ), shown in Fig. [16|, obtained with a = (x) 3 ^ 2 , 
(x) 2 , and e _5//3 (x) 2 , become reliable for the momentum transfer Q 2 > 95, 60, and 300 
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FIG. 16. Leading-twist NLO QCD results for Q 2 F 1V {Q 2 ) obtained with the 4>cz( x , Mf) dis- 
tribution amplitude and the choices of [i R given by Eqs. ( 4.19 ), with fip = Q 2 and (x) cz (Q 2 ) 



calculated according to ( [4.23 ). The solid curve (included for comparison) corresponds to the 
case = fip = Q 2 considered in the preceding subsection. The shaded area denotes the 
range of the total NLO prediction Q 2 F n (Q 2 ) where the upper limit corresponds to fi 2 R = ^"pms 
obtained from ( [5,11| ). 
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FIG. 17. The ratio F^\q 2 ) / F^\q 2 ) obtained with the 4>cz(x, fJ-p) distribution ampli- 
tude and the choices of fi R given by Eqs. ( |4.19| ), with fip = Q 2 and (x) cz {Q 2 ) calcu- 
lated according to ( |4.23 ). The solid curve (included for comparison) is for the result cor- 
responding to /J, R 



fip = Q 2 obtained in the preceding subsection, and the shaded area 
denotes the region of the predictions which correspond to as(^ 2 R ) <0.5 {ois{n 2 R ) <0.3), while 

\Fn\Q 2 )/FnHQ 2 ) 



< 0.3. 
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GeV 2 , respectively, if we demand that as(fi R ) < 0.3, and for Q 2 > 95, 35, and 50 GeV 2 , 
respectively, if as([i 2 R ) < 0.5 is regarded as a stringent enough requirement. 

Summarizing the above, one can say that contrary to the rather high value of Q 2 ~ 
500 GeV 2 required to obtain a reliable prediction assuming the (fi as (x, fi R ) distribution 
with fi R = Q 2 , one finds that by choosing the renormalization scale determined by the 
dynamics of the pion rescattering process, the size of the NLO corrections is significantly 
reduced and reliable predictions are obtained at considerably lower values of Q 2 , namely, 
for Q 2 < 100 GeV 2 . The same conclusion also applies to the results obtained with 
the 4>cz{x, fi 2 R ) distribution, for which choosing the renormalization scale related to the 
virtuality of the particles in the parton subprocess lowers the bound of the reliability of 
the results from Q 2 « 2400 GeV 2 to Q 2 < 100 GeV 2 . 



2. Theoretical uncertainty of the NLO predictions related to the renormalization scale 

ambiguity 

Unfortunately, at present we do not have at our disposal any absolutely reliable 
method of determining the "optimal" (or "correct") value of the renormalization scale 
for any particular order of PQCD. Our ignorance concerning the "optimal" value for fi R 
implies that any particular choice of this scale leads to an intrinsic theoretical uncertainty 



(error) of the perturbative results. Therefore, the NLO results displayed in Figs. [13] and 
|16| , being obtained with the four singled-out values of the renormalization scale [i R , 
contain theoretical uncertainty. In what follows we try to estimate this uncertainty. 
In order to do that, we have to make some assumptions, regarding the range of the 
renormalization scale ambiguity. 

As already mentioned in addition to the BLM method, two more renormalization 
scale-settings have been proposed: the FAC and PMS methods. All three of them are 
somewhat ad hoc and have no strong justification. Nevertheless, the principles underlying 
these methods are plausible, so that they give us at least a range of scale \i R which should 
be considered. 

Let UfaCi Upms, and /j,blm designate the scales determined by the BLM, PMS, and 
FAC scale setting methods, respectively. 

According to the FAC procedure, the scale fi R is determined by the requirement that 
the NLO coefficient in the perturbative expansion of F n (Q 2 ) vanishes, which, in our case, 
effectively reduces to solving the equation 

Fi 1 \Q 2 ^ 2 R = ^ 2 FAC )=0. (5.10) 

On the other hand, in the PMS procedure, one chooses the renormalization scale fx R 
at the stationary point of the truncated perturbative series for F n (Q 2 ). Operationally, 
this amounts to 



dF n (Q 2 ,y R ) 
d/4 



f^R — ^PMS ^ 



The BLM-determined scale is given by 



V 2 blm = e~ 5/3 (x) 2 Q 2 . (5.12) 
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The explicit expressions for F n (Q 2 ) and F^\Q 2 ) are given by (|5.3| - |5.7|) with fj/p — Q 



By solving Eqs. ( 5.10|) and (|5.11|) , and taking (|5.12|) into account we find that for the 
4>as{x, fip) distribution: 

fi 2 BLM = Q 2 m, (5.13a) 

/4ms ~Q 2 /18, (5.13b) 

» FAC ~Q 2 /1&, (5.13c) 

while, for the 4>cz(x,fi F ) distribution: 

Q 2 /79 <fi 2 BLM <Q 2 /77, (5.14a) 
Q 2 /59<fi 2 PMS <Q 2 /54, (5.14b) 
Q 2 /53 < fx 2 FAC < Q 2 /48 , (5.14c) 

for 50 GeV 2 < Q 2 < 100 GeV 2 . Therefore, we find that for both distributions 

I 1 BLM < I^PMS < ^FAC ■ 

(5.15) 

The FAC, PMS, and BLM scales, as evident from ( |5.13|) and (|5.14| ), are very close to 
each other, and the curves corresponding to the NLO prediction for Q 2 F n (Q 2 ) obtained 
with the FAC and PMS scales practically coincide with the dashed-dotted curves in Figs. 



13| and [T6| corresponding to the BLM scale. 

If the renormalization scale is interpreted as a "typical" scale of virtual momenta 
in the corresponding Feynman diagrams, then, despite the fact that we do not know 
what the "best" value of this scale is, it is (based on physical grounds and on the above 
considerations) reasonable to assume that it belongs to the interval ranging from Hblm 
to Q 2 . Namely, [i 2 BLM , being the lowest of the above considered scales and of the order 
of Q 2 /21 and Q 2 /80 for the (p as (x, fi 2 F ) and 4>cz(%, Hf) distributions, respectively, is a 
low enough to serve as the lower limit of the renormalization scale ambiguity interval. 
On the other hand, Q 2 is (too) high a scale and as such can safely be used for the upper 
limit of the same interval. 

In the following, then, instead of using any singled-out value, we vary the renormal- 
ization scale in the form 

li R = aQ 2 , (5.16a) 

where a is a continuous parameter 

ae[e~ 5/3 (x) 2 ,1]. (5.16b) 

Doing this will enable us to draw some qualitative conclusions concerning, first, the 
theoretical uncertainty related to the renormalization scale ambiguity, and, second, the 
effects that the inclusion of the NLO corrections has on the LO predictions. 

The LO result for the pion form factor is a monotonous function of the renormalization 
scale hr. Namely, all of the [Lr dependence of the LO prediction Q 2 Fl°\Q 2 ), as it is seen 
from (|5.3|) , is contained in the strong coupling constant o>s(l? r ). Thus, in accordance 
with Q2.2|), as [Ir decreases the LO result increases, and it increases without bound. 
In contrast to the LO, the NLO contribution Q 2 F^(Q 2 ), as evident from the explicit 
expression given in (|5.5| ), decreases (becomes more negative) with decreasing hr. Upon 
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adding up the LO and NLO contributions, we find that the full NLO result, as a function 
of stabilizes and reaches a maximum value for fi 2 R = fipMs- The values of the scale 
Upms are, for the 4> a s{x,Q 2 ) and (f>cz(x,Q 2 ) distributions given by ( |5.13b| ) and ( |5.14b[ ), 
respectively. 



If the renormalization scale continuously changes in the interval defined by ( |5.16| ), 
we find that the curves representing the LO and NLO results for Q 2 F W (Q 2 ) fill out the 



shaded regions in Figs. 12 and 13, for the <f) as (x,Q 2 ), and in Figs. |15] and [16] for the 
(pcz(x,Q 2 ) distribution. 

Next we turn to discuss the intrinsic theoretical uncertainty of the NLO prediction 
related to the renormalization scale ambiguity. Regarding this uncertainty, there is in 
fact no consensus on how to estimate it, or how to identify what the central value of 
Q 2 F n (Q 2 ) should be. 

Nevertheless, the simplest but still a good measure of this uncertainty is the quantity 
AF n (Q , /i m j n , H max ) — F n (Q , /i m j n ) — F^i^Q ,/i maa .), (5.17) 

i.e., the difference of the results for F n (Q 2 ) corresponding to the lower and the upper 
limit of the renormalization scale ambiguity interval ( |5.16| ). This quantity, therefore, for 
a given value of Q 2 , represents the "width" of the shaded regions in Figs. [12, O, [15], and 



IB] . In this sense, the shaded regions in these figures limited by the curves corresponding 
to fix = Q 2 and fi 2 R = fi 2 BLM essentially determine the theoretical accuracy allowed by 
the LO and NLO calculation. 

A glance at Figs. [12] and [13], displaying the LO and NLO predictions for Q 2 F 7T (Q 2 ) 
calculated with the 4> a s(%,Q 2 ) distribution, reveals that, compared to the LO, the NLO 
results exibit a much smaller renormalization scale dependence. The same holds true for 



the predictions depicted in Figs. [15| and |T^, based on the (ficz(%, Q 2 ) distribution. 



To be more quantitative, we thus find that, if, for the as (x, Q 2 ) distribution, at 
Q 2 = 50 GeV 2 (100 GeV 2 ), instead of /i 2 R = Q 2 one takes /i 2 R = fj? BLM = Q 2 /21, the LO 
result (Fig. [T2|) increases by 75% (64%), whereas the NLO result (Fig. |T3[) increases by 
20% (16%). Analogously, for the same values of Q 2 , but for the 4>cz(x,Q 2 ) distribution 
we find that taking fi R = n 2 BLM ~ Q 2 /79 instead [i 2 R = Q 2 the LO result increases by 
158% (125%), while the NLO result increases by 47% (34%). 

Therefore, the NLO corrections improve the situation because the terms in the NLO 
hard-scattering amplitude arise that cancel part of the scale dependence of the LO result. 

It should be pointed out that our estimate of the renormalization scale ambiguity 
interval given in ( |5.16| ) is very conservative, overestimating the theoretical uncertainty 
of the calculated NLO predictions. Namely, one could, almost at no risk replace Q 2 by 
) 3/2 Q 2 as the upper limit of the interval. If this is done, the dotted rather than solid 



curves would then provide the lower bound of the shaded regions in Figs. [L3| and [16 . 
Then, for Q 2 > 50 GeV 2 theoretical uncertainty of the NLO result for Q 2 F n (Q 2 ) turns 
out to be less than 5% for the (f) as {x, Q 2 ) and 8% for the 4>cz{x, Q 2 ) distribution. 

Before closing this subsection, a remark is appropriate. If the shaded areas in Figs. 



13| and [16] are displayed in the same figure they would not overlap. This implies that 
an unambiguous discrimination between the 4> as (x, Q 2 ) and <pcz(x, Q 2 ) distributions is 
possible, as soon as the data extending to higher values of Q 2 are obtained. 

Based on the above considerations, we may conclude that the inclusion of the NLO 
corrections stabilizes the LO prediction for the pion form factor by considerably reducing 
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the intrinsic theoretical uncertainty related to the renormalization scale ambiguity. This 
uncertainty for both distributions turns out to be of the order of a few percent. 



VI. SUMMARY AND CONCLUSIONS 



In this paper we have presented the results of a complete leading— twist NLO QCD 
analysis of the spacelike electromagnetic form factor of the pion at large momentum 
transfer. 

To clarify the discrepancies in the analytical expression for the hard— scattering am- 
plitude present in previous calculations, we have carefully recalculated the one— loop 
Feynman diagrams shown in Fig. |3|. Working in the MS renormalization scheme and 
employing the dimensional regularization method to treat all divergences (UV, IR, and 
collinear), we have obtained results which are in agreement with those of Refs. || (up 
to the typographical errors listed in ||) and 0. 

As nonperturbative input at the reference momentum scale of 0.5 GeV, we have used 
the four available pion distribution amplitudes defined by Eq. ([3.24J) and plotted in 
Fig. |[ the asymptotic distribution cj) as (x, fiF 2 ) and the three QCD sum— rule inspired 
distributions (ficz(x, Hf 2 ), 4>P2(x, Hf 2 ), and 4>ps(x, (if 2 )- The NLO evolution of these 
distributions has been determined using the formalism developed in Ref. ||20|| . 

By convoluting, according to Eq. ( |1 . 1|) , the hard-scattering amplitude with the pion 
distribution amplitude, both calculated in the NLO approximation, we have obtained the 
NLO numerical predictions for the pion form factor, for the four candidate distributions, 
and for several different choices of the renormalization and factorization scales, [ar and 
Lip- All the predictions have been obtained assuming n/ = 3 and Ajjg=0.2 GeV. 

We have first used the most simple choice of the scales where lir 2 = lif 2 = Q 2 ■ The 
results are summarized in Figs. ||] and |9] and Tables |II|-[V[ Although the predictions 
for the pion form factor Q 2 F n (Q 2 ) obtained with the four candidate distributions differ 
considerably, the lack of reliable experimental data, especially in the higher Q 2 region, 
does not allow us to confirm or discard any of them. The size of the NLO corrections 
and the size of the running coupling constant at given Q 2 can be used as indicators of the 
reliability of the perturbative treatment. Our numerical results based on the asymptotic 
distribution amplitude 4> as {x,LiF 2 ) differ from those of Ref. || (the difference is due to 
the different value of Ajjg). Thus, in contrast to Ref. U, where it was concluded that 
"reliable perturbative predictions can not be made until momentum transfers Q of about 
100 GeV are reached" , we have found that reliable predictions can already be made at 
momentum transfers of the order of 25 GeV. It has been shown that the inclusion of 
the (NLO) evolutional corrections only slightly influences the NLO prediction obtained 
assuming the <p as (x, li 2 f ) distribution. On the other hand, for the case of the end-point 
concentrated distributions, the evolutional corrections, both the LO and the NLO, are 
important. For the ii\ = fi 2 F = Q 2 choice of scales, the NLO corrections based on the 
cj) as (x, fMp), (j>cz(x, /Mp) 4>p2(x, fi F ), and 4>p3(x, fi F ) distributions, are large, which implies 
that one must demand that the momentum transfer Q be considerably larger than 10 
GeV before the corresponding results become reliable. 

In order to reduce the size of the NLO corrections and to examine the extent to 
which the NLO predictions for the pion form factor depend on the scales iar and lif, in 
addition to the simplest choice hr 2 = Hf 2 = Q 2 (which certainly is not best suited for 
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the process of interest), we have also considered the choices of /jlr and /if given by Eqs. 
( (4.1iJ| ) and ( |4.2(J| ) , respectively. The results appear to be insensitive to the choice of \ip 
considered. Using alternative choices for /x R , and modeling the pion with the (j) as (x, Q 2 ) 
and (f>cz(x,Q 2 ) distribution amplitudes, leads to the predictions shown in Figs. [13], [14], 
16] , and [T7|. The 4>p2(x, fiF 2 ) and 4>p3(x, /if 2 ) distributions are not separately considered, 



since the corresponding results are very similar to those obtained with (f>cz(x, (If 2 )- 

For a given distribution amplitude, the values of the pion form factor Q 2 F. 1T (Q 2 ) are 
rather insensitive to various choices of the scales This is evident from Figs. [13] and 



16], and is a reflection of the stabilizing effect that the inclusion of the NLO corrections 



has on the LO predictions. On the other hand, the ratio of the NLO corrections to 
the LO prediction, F« (Q 2 ) / 'F<® (Q 2 ) is very sensitive to the values of hr } as can be 
seen from Figs. 14 and |T7[ Requiring this ratio to be less than 0.3 and taking the less 
stringent condition on the value of the strong coupling a-s^p) < 0.5, we find that the 
predictions can be considered reliable for the momentum transfer Q < 10 GeV, provided 
the renormalization scale related to the average virtuality of the particles in the parton 
subprocess or given by the BLM scale, is used. 

Given the fact that we do not know what the "optimal" value of the renormalization 
scale is, choosing any particular value for this scale introduces a theoretical uncertainty in 
the NLO predictions. Based on a reasonable guess of the renormalization scale ambiguity 
interval, we have estimated this uncertainty to be less than 10%. 

The difference between the absolute predictions based on the (fi as (x, fi 2 F ), 4>cz(x,fj, F ) 
(4>P2(x, ii F ))i an d (pP3(x,fi F ) distributions is large enough to allow an unambiguous ex- 
perimental discrimination between them, as soon as the data extending to higher values 
of Q become available. 

In conclusion, the results of the complete leading-twist NLO QCD analysis, which 
has been carried out in this paper, show that reliable perturbative predictions for the 
pion electromagnetic form factor with all the four distribution amplitudes considered 
can already be made at a momentum transfer Q < 10 GeV. The theoretical uncertainty 
related to the renormalization scale ambiguity, which constitutes a reasonable range of 
physical values, has been shown to be less than 10%. To check our predictions and to 
choose between the distribution amplitudes, it is necessary that experimental data at 
higher values of Q 2 are obtained. 



ACKNOWLEDGMENTS 



The authors would like to thank A. V. Radyushkin for pointing out an error present 
in the original version of the manuscript, and P. Kroll for useful suggestions. This work 
was supported by the Ministry of Science and Technology of the Republic of Croatia 
under Contract No. 00980102. 



37 



REFERENCES 

[1] S. J. Brodsky and G. P. Lepage, Phys. Lett. 87B, 359 (1979); Phys. Rev. Lett. 43, 
545 (1979); Phys. Rev. Lett. 43, 1625 (1979) (E); G. P. Lepage and S. J. Brodsky, 
Phys. Rev. D 22, 2157 (1980). 

[2] A. V. Efremov and A. V. Radyushkin, Theor. Mat. Phys. 42, 97 (1980); Phys. Lett. 
94B, 245 (1980). 

[3] A. Duncan and A. H. Mueller, Phys. Lett. 90B, 159 (1980); Phys. Rev. D 21, 1636 
(1980). 

[4] S. J. Brodsky and G. P. Lepage, in Perturbative QCD, edited by A. H. Mueller (World 
Scientific Publishing Co., Singapore, 1989); V. L. Chernyak, in High pr Physics and 
Higher Twists, Proceedings of the Conference, Paris, France, 1988, edited by M. 
Benayoun, M. Fontannaz, and J. L. Narjoux, [Nucl. Phys. Proc. Suppl. 7B, 297 
(1989)]. 

[5] R. D. Field, R. Gupta, S. Otto, and L. Chang, Nucl. Phys. B186, 429 (1981). 

[6] F.-M. Dittes and A. V. Radyushkin, Yad. Fiz. 34, 529 (1981), [Sov. J. Nucl. Phys. 

34, 293 (1981) ]. 
[7] M. H. Sarmadi, Ph. D. thesis, University of Pittsburgh, 1982. 

[8] R. S. Khalmuradov and A. V. Radyushkin, Yad. Fiz. 42, 458 (1985), [Sov. J. Nucl. 

Phys. 42, 289 (1985)]. 
[9] E. Braaten and S.-M. Tse, Phys. Rev. D 35, 2255 (1987). 
[10] E. P. Kadantseva, S. V. Mikhailov, and A. V. Radyushkin, Yad. Fiz. 44, 507 (1986), 

[Sov. J. Nucl. Phys. 44, 326 (1986) ]. 
[11] F. D. Aguila and M. K. Chase, Nucl. Phys. B193, 517 (1981). 
[12] E. Braaten, Phys. Rev. D 28, 524 (1983). 
[13] B. Nizic, Phys. Rev. D 35, 80 (1987). 

[14] N. Isgur and C. H. Llewellyn Smith, Phys. Rev. Lett. 52, 1080 (1984); Phys. Lett. 

217B, 535 (1989); Nucl. Phys. B 317, 526 (1989). 
[15] A. V. Radyushkin, Nucl. Phys. A 532, 141 (1991). 

[16] H.-N. Li and G. Sterman, Nucl. Phys. B 381, 129 (1992); D. Tung and H.-N. Li, 

Chin. J. Phys. 35, 651 (1997). 
[17] R. Jakob and P. Kroll, Phys. Lett. 315B, 463 (1993); ibid. 319B, 545 (E) (1993). 
[18] F.-M. Dittes and A. V. Radyushkin, Phys. Lett. 134B, 359 (1984); M. H. Sarmadi, 

ibid. 143B, 471 (1984); G. R. Katz, Phys. Rev. D 31, 652 (1985); S. V. Mikhailov 

and A. V. Radyushkin, Nucl. Phys. B254, 89 (1985). 
[19] S. V. Mikhailov and A. V. Radyushkin, Nucl. Phys. B273, 297 (1986). 
[20] D. Miiller, Phys. Rev. D 49, 2525 (1994); Phys. Rev. D 51, 3855 (1995). 
[21] G. Curci, W. Furmanski, and R. Petronzio, Nucl. Phys. B175, 27 (1980). 
[22] S. J. Brodsky, P. Damgaard, Y. Frishman, and G. P. Lepage, Phys. Rev. D 33, 1881 

(1986). 

[23] V. L. Chernyak and A. R. Zhitnitsky, Phys. Rep. 112, 173 (1984). 
[24] G. Martinelli and C. Sacharadja, Phys. Lett. 190B, 151 (1987); Phys. Lett. 217B, 
319 (1989). 

[25] E. G. Floratos, D. A. Ross, and C. T. Sachrajda, Nucl. Phys. B129, 66 (1977); ibid. 

B 139, 545 (E) (1978). 
[26] A. Gonzales-Arroyo, C. Lopez, and F. J. Yndurain, Nucl. Phys. B153, 161 (1979). 
[27] G. R. Farrar, K. Huleihel, and H. Zhang, Nucl. Phys. B349, 655 (1991). 
[28] S. J. Brodsky, G. P. Lepage, and P. B. Mackenzie, Phys. Rev. D 28, 228 (1983). 



38 



[29] G. Grunberg, Phys. Lett. 95B, 70 (1980); Phys. Lett. HOB, 501 (1982); Phys. Rev. 
D 29, 2315 (1984). 

[30] P. M. Stevenson, Phys. Lett. 100B, 61 (1981); Phys. Rev. D 23, 2916 (1981); Nucl. 

Phys. B203, 472 (1982); Nucl. Phys. B231, 65 (1984). 
[31] G. Kramer and B. Lampe, Z. Phys. A 339, 189 (1991). 

[32] G. Parisi and R. Petronzio, Phys. Lett. 94B, 51 (1980); A. C. Mattingly and P. M. 
Stevenson, Phys. Rev. D 49, 437 (1994); V. N. Gribov, Lund Report No. LU-TP 
91-7, 1991 (unpublished); K. D. Born, E. Laermann, R. Sommer, P. M. Zerwas, and 
T. F. Walsh, Phys. Lett. 329B, 325 (1994); J. M. Cornwall, Phys. Rev. D 26, 1453 
(1982); A. Donnachie and P. V. Landshoff, Nucl. Phys. B311, 509 (1989). 

[33] S. Peris and E. de Rafael, Nucl. Phys. B500, 325 (1997). 

[34] M. B. Gay Ducati, F. Halzen, and A. A. Natale, Phys. Rev. D 48, 2324 (1993). 
[35] J. Bebek et al, Phys. Rev. D 17, 1693 (1978). 
[36] S. R. Amendolia et al., Nucl. Phys. B277, 168 (1986). 
[37] C. E. Carlson and J. Milana, Phys. Rev. Lett. 65, 1717 (1990). 
[38] S. Dubnicka and L. Martinovic, Phys. Rev. D 39, 2079 (1989); J. Phys. G 15, 1349 
(1989). 

[39] A. V. Radyushkin, Acta Phys. Polon. B26, 2067 (1995); S. Ong, Phys. Rev. D 52, 
3111 (1995); P. Kroll and M. Raulfs, Phys. Lett. 387B, 848 (1996). 



39 



